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Using a sample of 9.6 x 10 6 BB meson pairs collected with the CLEO de- 
' tector, we have fully reconstructed 135 B° ^>W Kg and 526 5+ -► V W ^ + 

candidates with very low background. We fitted the K invariant mass dis- 
tributions of these B meson candidates and measured the masses of the neutral 
and charged B mesons to be M(B°) = 5279.1±0.7[stat]±0.3[syst] MeV/c 2 and 
M(B+) = 5279.1 ±0.4[stat] ±0.4[syst] MeV/c 2 . The precision is a significant 
' improvement over previous measurements. 
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The previous measurements of the B meson masses at e + e colliders operating at T(45) 
energy Jl| were obtained from fits to the distributions of the beam-constrained B mass, 

defined as M bc = \J E^ cam — p 2 (B), where p(B) is the absolute value of the B candidate 
momentum. Substitution of the beam energy for the measured energy of the B meson 
candidate results in a significant improvement of the mass resolution, therefore the beam- 
constrained mass method is the technique of choice for the M(B°) — M(B + ) mass difference 
measurement. However, the precision of the measurement of the absolute B° and B + meson 
masses is limited by the systematic uncertainties in the absolute beam energy scale and in 
the correction for initial state radiation. For this measurement we selected B° — > ip^ Kg 
and B + — > tp^ K + candidates f2j, reconstructing ip^ — > £ + £~ and K s — > decays. We 

used both e + e~ and modes for the ip^ reconstruction. We then determined the B° 

and B + meson masses by fitting the ip^ K s and K + invariant mass distributions. The 
main reasons to use the ip^K rather than more copious D^nn final states are, first, that the 
background is very low; second, that the J/ip and ip{2S) mesons are heavy, and their masses 
are very well measured 0]. As discussed below, constraining the reconstructed J/ip and 
ip(2S) masses to their world average values makes our B mass measurement insensitive to 
imperfections in the lepton momentum reconstruction. By comparing the beam-constrained 
B mass to the B° and B+ mass values obtained in our measurement, one could set the 
absolute beam energy scale at the e + e~ colliders operating in the T(4S) energy region. 

The data were collected at the Cornell Electron Storage Ring (CESR) with two configu- 
rations of the CLEO detector, called CLEO II @ and CLEO II. V §. The components of 
the CLEO detector most relevant to this analysis are the charged particle tracking system, 
the Csl electromagnetic calorimeter, and the muon chambers. In CLEO II, the momenta 
of charged particles are measured in a tracking system consisting of a 6-layer straw tube 
chamber, 10-layer precision drift chamber, and 51-layer main drift chamber, all operating 
inside a 1.5 T solenoidal magnet. The main drift chamber also provides a measurement of the 
specific ionization, dE/dx, used for particle identification. For CLEO II. V, the straw tube 
chamber was replaced with a 3-layer silicon vertex detector. The muon chambers consist of 
proportional counters placed at various depths in the steel absorber. Track fitting was per- 
formed using a Kalman filtering technique, first applied to track fitting by P. Billoir |J. The 
track fit sequentially adds the measurements provided by the tracking system to correctly 
take into account multiple scattering and energy loss of a particle in the detector material. 
For each physical track, separate fits are performed using different particle hypotheses. 

For this measurement we used 9.1 fb _1 of e + e~ data taken at the T(4S") energy and 
4.4 fb _1 recorded 60 MeV below the T(4S') energy. Two thirds of the data used were 
collected with the CLEO II. V detector. All of the simulated event samples used in this 
analysis were generated with a GEANT-based simulation of the CLEO detector response 
and were processed in a similar manner as the data. 

Electron candidates were identified based on the ratio of the track momentum to the 
associated shower energy in the Csl calorimeter and specific ionization in the drift chamber. 
We recovered some of the bremsstrahlung photons by selecting the photon shower with 
the smallest opening angle with respect to the direction of the track evaluated at the 
interaction point, and then requiring this opening angle to be smaller than 5°. For the 
— > reconstruction one of the muon candidates was required to penetrate the steel 
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absorber to a depth greater than 3 nuclear interaction lengths. We relaxed the absorber 
penetration requirement for the second muon candidate if it pointed to the end-cap muon 
chambers and its energy was too low to reach a counter. For these muon candidates we 
required the ionization signature in the Csl calorimeter to be consistent with that of a 
muon. 

We extensively used normalized variables, taking advantage of well-understood track 
and photon-shower four-momentum covariance matrices to calculate the expected resolution 
for each combination. The use of normalized variables allows uniform candidate selection 
criteria to be applied to the data collected with the CLEO II and CLEO II. V detector 
configurations. The e + (7)e~(7) and /U + /i~ invariant mass distributions for the ipw — ► £ + £~ 
candidates in data are shown in Fig. [I]. We selected the ip('> — > £ + £~ signal candidates 
requiring the absolute value of the normalized invariant mass to be less than 3. The average 
£ + £~ invariant mass resolution is approximately 12 MeV/c 2 . For each if)W candidate we 
performed a fit constraining its mass to the world average value |3|]. This mass-constraint 
fit improves the ip^ energy resolution almost by a factor of 4 and the absolute momentum 
resolution by 30%. 
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FIG. 1. (a) ift^ — > e + e~ and (b) ip^ — > [i^ [i~ candidates in data. The solid line represents the 
data taken at T(45) energy; the dashed line represents the scaled off resonance data showing the 
level of background from non-Si? events. The shaded parts of the histograms represent the 
candidates with the absolute value of the normalized invariant mass less than 3. 

The Kg candidates were selected from pairs of tracks forming well-measured displaced 
vertices. The daughter pion tracks were re-fitted taking into account the position of the dis- 
placed vertex and were constrained to originate from the same spatial point. The resolution 
in 7T + 7r _ invariant mass is approximately 4 MeV/c 2 . After requiring the absolute value of 
the normalized tt + it~ invariant mass to be less than 3, we performed a fit constraining the 
mass of each K s candidate to the world average value ||. 
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The B — > ip^ K candidates were selected by means of two observables. The first ob- 
servable is the beam-constrained B mass M^ c = \J E^ CSLm — p 2 (B). The resolution in M bc 
for the B — ► ip® K candidates is approximately 2.7 MeV/c 2 and is dominated by the 
beam energy spread. We required \M^ C — 5280 MeV/c 2 |/<r(Mb c ) < 3. The requirement 
on Mbc equivalent to a requirement on the absolute value of the B candidate momentum, 
is used only for background suppression, and does not bias the B mass measurement. The 
second observable is the invariant mass of the ip^ K system. The average resolutions in 
M(> (/) Kfj and M(^ (/) K + ) are respectively 8 MeV/c 2 and 11 MeV/c 2 . The M(^ (/) K) dis- 
tributions for the candidates passing the beam-constrained B mass requirement are shown 
in Fig. |. To select signal candidates, we required |M(^ {/) Kg) - 5280 MeV/c 2 |/a(M) < 4 
and |M(> (/) K + ) - 5280 MeV/c 2 |/cr(M) < 3; the allowed invariant mass intervals are suffi- 
ciently wide not to introduce bias in the B mass measurement. These selections yielded 135 
B° -> V W K s candidates: 125 in the B° -> J/^lTg mode and 10 in the B° -> ^(2S r ) ATS} 
mode. We estimated the background to be 0.13i[J;o5 events. The selections yielded 526 
B+ _> ^0^+ candidates: 468 in 5+ -> J/^K + mode and 58 in B+ -> ^(25)^+ mode. 
The background from 5 + — ► t/>W 7r + decays was estimated to be 0.9 ± 0.3 events, whereas 
all other background sources were estimated to contribute 2. 3i ;s events. The backgrounds 
were evaluated with simulated events and the data recorded at the energy below the BB pro- 
duction threshold. We discuss the systematics associated with background below, together 
with other systematic uncertainties. 
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FIG. 2. The invariant mass distributions for (a) B° -> V W and (b) B + -> V> (/) ^ + candidates 
passing the beam-constrained -B mass requirement. The shaded parts of the histograms represent 
the candidates selected for the B mass measurement fits by requiring that the absolute value of the 
normalized invariant mass to be less than 4 in (a) and 3 in (b). 

The B meson masses were extracted from the ip^ K s and K + invariant mass distri- 
butions with an unbinned likelihood fit. The likelihood function is 
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where Mj is the invariant mass of a ip^ K combination, Oi is the calculated invariant mass un- 
certainty for that tp^ K combination, and G(x\a) = l/(y/2n(r) exp(— x 2 /2a 2 ). The product 
is over the B° or B + meson candidates. The parameters of the fit are the B meson mass M(B) 
and a global scale factor S that modifies the calculated invariant-mass uncertainties (Tj. The 
integration limits of the normalization integral in the denominator correspond to the signal 
regions defined for the ip^ K invariant mass distributions: [5280 MeV/c 2 ±4<j(M)] for B° and 
[5280 MeV/c 2 ±3(x(M)] for B + candidates. From the fits to the tp^ K% invariant-mass distri- 
bution, we obtained M(B°) = 5278.97±0.67 MeV/c 2 , S = 1.24±0.08, and the correlation co- 
efficient p(M(B°), S) = -0.013. For ^ (/) K + we obtained M(B + ) = 5279.50 ± 0.41 MeV/c 2 , 
S = 1.09 ± 0.04, and p(M(B + ), S) = —0.015. The values of the scale factor S and uncer- 
tainty in M(B) returned by the fits are in good agreement with the values obtained from 
simulated events. 

Table I lists the bias corrections together with associated systematic uncertainties, which 
we will discuss below. 

TABLE I. Bias corrections and associated systematic uncertainties. The total systematic un- 
certainty was obtained by adding in quadrature all the uncertainties listed in the table. 



Source of bias Correction and uncertainty 



(MeV/c 2 ) 

M (B°) M{B + ) 

(i) Bias observed for simulated events +0.08 ± 0.08 -0.13 ±0.13 

(ii) Background ±0.05 ±0.17 

(iii) B mass likelihood fit ±0.14 ±0.09 

(iv) ipW four-momentum measurement ±0.05 ±0.06 

(v) Kg four-momentum measurement ±0.28 

(vi) K + momentum measurement —0.32 ± 0.34 

(vii) Detector misalignment negligible negligible 



Total correction and systematic uncertainty ±0.08 ± 0.33 —0.45 ± 0.42 



(i) Measuring B masses using simulated events. - - Applying the same procedure as in 
the data analysis, we measured the B° and B + mass using 30286 B° -> ^ (,) K° s and 34519 
B + — > K + candidates reconstructed from a sample of simulated events. We obtained 
M(B°) - MjjT* = -0.08 ± 0.04 MeV/c 2 and M(B + ) - M^T* = ±0.13 ± 0.05 MeV/c 2 . We 
applied ±0.08 MeV/c 2 and —0.13 MeV/c 2 corrections to the B° and B + mass values and 
assigned 100% of those corrections as systematic uncertainties. 

(ii) Background. - - The estimated mean background for the B° candidates is 0. 131^05 
events, we therefore conservatively assumed the probability of finding a background event 
in our sample to be 22%. The B° background candidates are expected to be uniformly 
distributed across the B° mass signal region. We performed the B° mass fits excluding the 
one candidate with the highest or the lowest normalized Kg invariant mass; the largest 
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observed B° mass shift was 0.25 MeV/c 2 . We multiplied this 0.25 MeV/c 2 shift by the 22% 
probability of having a background event in our sample and assigned 0.05 MeV/c 2 as the 
systematic uncertainty in B° mass due to background. For the B + signal, the background 
from B + — > if>^'> 7r + decays was estimated to be 0.9 ± 0.3 events, all other background 
sources were estimated to contribute 2.3^5 events. The B + background candidates, with 
the exception of B + — > if>('> n + events, are expected to be uniformly distributed across the 
B mass signal region. The B + — > if>('' ix + events reconstructed as B + — > if)^ K + produce 
high tp('> K + invariant mass. We performed the B + mass fits excluding 4 candidates with 
the highest or the lowest normalized if>^ K + invariant mass and assigned the largest shift 
of the measured B + mass (0.17 MeV/c 2 ) as the systematic uncertainty in B + mass due to 
background. 

(iii) B mass likelihood fit. - - We studied the systematics associated with the unbinned 
likelihood fit procedure by changing the fit function from a Gaussian to sum of two Gaussians. 
We also allowed the fit to determine different scale factors 5* for the candidates coming from 
CLEO II and CLEO II. V data, or for the candidates with if>^ — > e + e~ and if>('> — > 

We assigned the largest shift of the measured B mass as the systematic uncertainty. 

(iv) if>('> four-momentum measurement. — Even if the B mass measurements using the 
simulated events show negligible bias, a bias in the measurement is still in principle possible 
because of the uncertainty in the absolute magnetic field scale, an imperfect description of 
the detector material used by the Billoir fitter, or detector misalignment. For the if)*' four- 
momentum measurement, these systematic effects along with the systematics associated with 
bremsstrahlung are rendered negligible by the the ipW mass-constraint fit. The measured 
position of the J /if) mass peak allows a reliable evaluation of the possible bias in the lepton 
momentum measurement. We measured the positions of the J /if) — > fi + fi~ and J/if) —>■ e + e~ 
peaks by fitting the inclusive and e + ( r y)e~( r y) invariant mass distributions. In these 
fits we used the signal shapes derived from a high-statistics sample of simulated J /if) — > £ + £~ 
events generated with the J /if) mass of 3096.88 MeV/c 2 |J. In simulated J /if) — > events, 
the reconstruction procedure introduces a bias of less than 0.03 MeV/c 2 in the measured 
J /if) mass. We found that the J /if) — » peak was shifted by +0.5 ± 0.2 MeV/c 2 in 
data compared to simulated events; the corresponding value of the J /if) — > e + e~ peak shift 
was +0.7 ± 0.2 MeV/c 2 . A +0.5 MeV/c 2 shift corresponds to overestimation of the lepton 
absolute momenta by approximately 0.02%. A variation of the lepton absolute momenta 
by 0.1% produced a shift of less than 0.02 MeV/c 2 in the measured B mass. We therefore 
neglected the systematic uncertainty associated with the lepton momentum measurement. 
In addition, we varied the world average J /if) and if)(2S) mass values used in the mass- 
constraint fits by one standard deviation ||; the resulting 0.05 MeV/c 2 and 0.06 MeV/c 2 
shifts of the measured B° and B + masses were assigned as systematic uncertainties. 

(v) Kg four-momentum measurement — The systematic uncertainty of our B mass mea- 
surement is dominated by a possible bias in the kaon four- momentum measurement. The 
measured position of the Kg mass peak allows a reliable evaluation of the possible bias in the 
Kg four- momentum measurement. We selected inclusive K s candidates satisfying the same 
K s selection criteria as in the B° — > ipw K s analysis; the momenta of the selected inclusive 
Kg candidates were further restricted to be from 1.55 to 1.85 GeV/c, which corresponds 
to the momentum range of the K s mesons from B° — > J/ipK® decays. Using this sample, 
we measured the mean reconstructed K° mass to be within 10 keV/c 2 of the world average 
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value of 497.672 ± 0.031 MeV/c 2 ||. However, we also observed a ±40 keV/c 2 variation of 
the measured mean K° mass depending on the radial position of the K s decay vertex. To 
assign the systematic uncertainty, we conservatively took the K° mass shift to be 40 keV/c 2 
and added in quadrature the 30 keV/c 2 uncertainty in the world average K° mass to obtain 
a total shift of 50 keV/c 2 . This 50 keV/c 2 variation in the measured Kg mass could be 
obtained by varying each daughter pion's momentum by 0.018%; the resulting variation of 
the measured B° mass was 0.26 MeV/c 2 , which we assigned as a systematic uncertainty due 
to the K s four-momentum measurement. This uncertainty in M(B°) has a contribution 
from the uncertainty in the world average value of the K° mass, which partially limited the 
precision of our Kg mass peak position measurement. In addition, we varied by one standard 
deviation the world average K° mass value used for the K s mass-constraint fit; the resulting 
0.04 MeV/c 2 variation of the measured B° mass was added to the systematic uncertainty. 

(vi) K + momentum measurement. — Comparing the momentum spectra of the muons 
from inclusive J/ip decays and the kaons from B + — > tp^ K + decays, we concluded that 
J/ij) — > n + ji~ decays provide excellent calibration sample for the study of the systematic 
uncertainty associated with the K + momentum measurement. As discussed above, the ob- 
served +0.5 ± 0.2 MeV/c 2 shift of the J/ip — > mass peak corresponds to a systematic 
overestimation of the muon momenta by 0.02%. We decreased the measured K + momenta 
by 0.02%, which resulted in a -0.32 MeV/c 2 shift of the measured B + mass. We applied 
a —0.32 MeV/c 2 correction to our final result and assigned 100% of the correction value as 
the systematic uncertainty. The ionization energy loss for muons from inclusive J/^'s dif- 
fers slightly from that for kaons from B + — > ipW K + decays. To account for the systematic 
uncertainty due to this difference, we measured the B + mass using the pion Billoir fit hy- 
pothesis for kaon tracks. The resulting shift (0.08 MeV/c 2 ) was added in quadrature to the 
systematic uncertainty. Because of acceptance of the muon chambers, the muons pointing 
to the end-cap region of the detector are under-represented in comparison with the kaons 
from the B + — > if)W K + decays. The tracks with low transverse momentum are more likely 
to be affected by the magnetic field inhomogeneity, thus providing an additional source of 
systematic bias, which will not be taken into account by studying J/tp — > fi + fi~ decays. 
However, if a K + track has low transverse momentum, then its track parameters are poorly 
measured, and the mass fit naturally assigns a low weight to this B + candidate. We studied 
the possible systematic bias both by varying the measured K + momentum by 0.1% for the 
K + tracks with | cos^l > 0.8, where 9 is the angle between a track and the beam direction, 
and by excluding these low angle tracks altogether. The largest shift (0.08 MeV/c 2 ) was 
added in quadrature to the systematic uncertainty. 

(vii) Detector misalignment. — The detector misalignment effects were studied with high- 
momentum muon tracks from e + e~ — > fi + fi~ events. We measured the mean of the transverse 
momentum difference between the /i + and \i~ tracks. We also studied the dependence of the 
sum of the fi + and fi~ momenta on azimuthal angle and polar angle 9 of the fi + track. 
We parametrized our findings in terms of an average as well as <fi- and ^-dependent false 
curvature. We varied the measured curvature of the signal candidate tracks according to 
these parametrizations and found the detector misalignment effects to be negligible for our 
B mass measurements. 

In conclusion, we have determined the masses of neutral and charged B mesons with 
significantly better precision than any previously published result |§. We obtained M(B°) = 
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5279.1 ±0.7[stat]±0.3[syst] MeV/c 2 and M(B+) = 5279.1 ±0.4[stat] ±0.4[syst] MeV/c 2 . The 
systematic uncertainties for the M(B°) and M(B + ) measurements are independent except 
for the small common uncertainties due to the imperfect knowledge of the J /if) and if)(2S) 
masses (item (iv) in Table I). Combining our M(B°) and M(B + ) measurements with the 
world average value of the mass difference M(B°) - M{B+) = 0.34 ± 0.32 MeV/c 2 [§, we 
obtained M(B°) = 5279.2 ± 0.5 MeV/c 2 and M(B+) = 5278.9 ± 0.5 MeV/c 2 . Although 
these M(B°) and M(B + ) values are more precise than the results given above, obviously 
they are strongly correlated: the correlation coefficient is p(M(B°), M(B + )) = 0.81. 
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